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ABSTRACT 

Solar prominences are an important tool for studying the structure and evolution of the coronal magnetic field. 
Here we consider so-called "hedgerow" prominences, which consist of thin vertical threads. We explore the 
possibility that such prominences are supported by tangled magnetic fields. A variety of different approaches 
are used. First, the dynamics of plasma within a tangled field is considered. We find that the contorted shape 
of the flux tubes significantly reduces the flow velocity compared to the supersonic free fall that would occur 
in a straight vertical tube. Second, linear force-free models of tangled fields are developed, and the elastic 
response of such fields to gravitational forces is considered. We demonstrate that the prominence plasma can 
be supported by the magnetic pressure of a tangled field that pervades not only the observed dense threads but 
also their local surroundings. Tangled fields with field strengths of about 10 G are able to support prominence 
threads with observed hydrogen density of the order of lO" cm"^. Finally, we suggest that the observed vertical 
threads are the result of Rayleigh-Taylor instability. Simulations of the density distribution within a prominence 
thread indicate that the peak density is much larger than the average density. We conclude that tangled fields 
provide a viable mechanism for magnetic support of hedgerow prominences. 
Subject headings: MHD — Sun: corona — Sun: magnetic fields — Sun: prominences 



1. INTRODUCTION 

Solar prominences (a.k.a. filaments) are relatively cool 
structures embedded in the million-degree corona at heights 
well above the chromosphere (see reviews by Hiravama 1985; 
Zirkdll989t lPriestlll990t lTandber g-Hansseni il995c iHeinzell 
Above the solar limb, prominences appear as bright 
structures against the dark background, but when viewed as 
filaments on the solar disk they can be brighter or darker than 
their surroundings, depending on the bandpass used to ob- 
serve them. Magnetic fields are thought to play an impor- 
tant role in supporting the prominence plasma against gravity, 
and in insulating it from the surrounding hot corona. Most 
quiescent prominences exhibit intricate filamentary struc- 
tures that evolve continually due to plas ma flows and heat- 
ing/cooling processes (see examples in Menzel & Wolbach 
1960; Engvold 1976; Malherbe 1989; Leroy 1989; Martin 
19981) . In some cases the threads appear to be mostly hori- 
zontal, while in other cases they are clearly radially oriented 
("hedgerow" prominences). Figure[T]shows several examples 
of prominences observed in Ha at the Big Bear Solar Obser- 
vatory (BBSO) and the Dutch Open Telescope (DOT). The ex- 
amples in Figs.[T^ and[TJ) show mainly vertical threads, while 
the prominence in Fig.[T]; shows horizontal threads. Off-limb 
observations in He II 304 A indicate that there are higher alti- 
tude parts that are optically thin in Ha and therefore not visi- 
ble on the disk (or at least have not been clearly identified in 
disk observations). Figure [TJl shows a prominence above the 
solar limb a s observed in He II 3 04 A with the SECCHI/EUVI 
instrument (iHoward et al.ll20(j7h on the STEREO spacecraft. 
The upper parts of the prominence consist of vertical threads 
with an intricate fine-scale structure. Movie sequences of qui- 
escent and erupting prominences can be found at the STEREO 
website'. 

Prominence plasma is highly dynamic, ex- 
hibiting horizontal and vertical motions of or- 

' http://stereo.gsfc.nasa.gov/gallery/selects.shtml 




Fig. 1. — Quiescent prominences above the solar limb: (a) Ha image of 
prominence observed at Big Bear Solar Observatory (BBSO), November 22, 
1995; (b) Ha, BBSO, 1970; (c) Ha, Dutch Open Telescope, September 15, 
2006; (d) He II 304 A, STEREO/EUVl, 2008 April 20 at 00:06 UT. 



der 10 -70 km s:^ (Menzel & Wolbach 19601 

Engvoldl 119761; IZirker. Engvold & Martin 1991 
Kucera. Tovar & De Pointieul 2003; Lin, Engvold & Wiikl 
'2003^. lOkamoto et al.1 120071 llergeret al.i i2008t iChae et al.1 
2008). Recent high-resolution observations of filaments on 
the solar disk indicate that they consist of a collection of very 
thin threads with widths of about 200 km, at the limit of reso- 
lution of present-da y telescopes dOn . Engvold & Wiik 20031 
'Lin et al.l l2005aS iLin. Martin & Engvoldi i2008t iLin et all 
2008). Individual threads have lifetimes of only a few min- 
utes, but the filament as a whole can live for many days. It 
seems likely that these thin threads are aligned with the local 
magnetic field. High-resolution images of prominences above 
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the Hmb have been obtained with the S olar Optical Telescope 
(SOT) onboard Hinode. For example. lOkamoto et alJ (l2007h 
observed horizontal threads in a prominence near an active 
region, and st u died the oscillatory motions of these threads. 
iHeinzel et alJ (2008) observed a hedgerow prominence 
consisting of many thin vertical threads, and they used 
multi-wavelength observations to estimate the amount of 
absorption and "emis sivity blocking" in t he prominence and 
surrounding cavity. iBerger et al.l (l2008h observed another 
hedgerow prominence and found that the prominence sheet 
is structured by both bright quasi- vertical threads and dark 
inclusions. The bright structures are downflow streams with 
velocity of about 10 km s"', and the dark inclusions are 
highly dynamic upflows with velocity of about 20 km s"'. 
The downflow velocities are much less than the free-fall 
speed, indicating that the plasma is so mehow being supported 
against gravity. iBerger et al.l (l2008l) proposed that the dark 
plumes conta in relatively hot p lasma that is driven upward 
by buoyancy. IChae et alT(l200^ ) observed a persistent flow of 
Ha emitting plasma into a prominence from one side, leading 
to the formation of vertical threads. They suggested that 
the vertical threads are stacks of plasma supported against 
gravity by the sagging of initially horizontal magnetic field 
lines. 

Direct measurements of the prominence magnetic field 
can be o btained using spec tro-polarimetr y (see reviews 
bv iLerovl 1989: Paletou & Aulanied 120031: iPaletoul 120081: 
iLopez Ariste & Aulanier .20071) . A comprehensive effort to 
measure prominence magnetic fields was conducted in the 
1970's and early 1980's using the facilities at Pic du Midi 
(France) and Sacramento Peak Observatory (USA). This work 
showed that (1) the magnetic field in quiescent prominences 
has a strength of 3-15 G; (2) the field is mostly horizon- 
tal and makes an angle of about 40° with respect to the 
long axis of the prominence (Lerovl l 19891: iBommier & Lerovl 
[1998; Paletou & Aulanier 2003); (3) the field strength in- 
creases slightly with height, indicating the presence of dipped 
field lines; (4) most prominences have inverse polarity, i.e., 
the component of magnetic field perpendicular to the promi- 
nence axis has a direction opposite to that of the potential 
field. These earlier data likely included a variety of quiescent 
prominences, some with predominantly horizontal threads, 
others with more v ertical threads. In more recent work, 
iPaletou et al.l (1200 ll) reported full-Stokes observations of a 
limb prominence in He I 5876 A (He D3), and derived mag- 
netic field strengths of 30-45 G, s omewhat lar ger than the 
values reported in earlier studies. ICasini et"an 42003) pub- 
lished the first vector-field map of a prominence with a spatial 
resolution of a few arcseconds. They found that the average 
magnetic field in this prominence is mostly horizontal with a 
strength of about 20 G and with the magnetic vector pointing 
20° to 30° off the prominence axis, consistent with the ear- 
lier studies. However, the map also shows clearly organized 
patches where the magnetic field i s significantly stronger 
than average, up to 80 G (also see Wiehr & Bianda 2003 
Lopez Ariste & Casini 2002, 2003; Lopez Ariste & Aulanier 
20071: ICasini. Bevilacqua & Lopez Aristell2005l) . It is unclear 



ported against gravity? Many authors have suggested 
that quiescent prominences are embedded in large-scale 
flux ropes that lie horizontally above the polarity in- 
version line d Kuperus & Raadu 1974; Pneuman 1983t 
van Ballegooiien & Martens 1989; Priest, Hood & Anzej 
1989; Rust &Kumar . 1_994; .Lqw _& Hund hausen 19951: 



Aulan ier et al.1 ll 9981: IChae et al"n200ir " lGibIon & Fanl 
i2006i: iDudik et al.l " l20m The prominence plasma 
is thought to be located near the dips of the helical 
field lines. The magnetic field near the dips may be 
deformed by the wei gh t of th e prominence plasma 
(iKippenhahn & Schliited \m3k |Lok&£etrid I200I 
iPetrie & Low! 120051: iHeinzel. Anzer & Gunad l2005l) . Oth- 



how these patches are related to the fine threads seen at higher 
spatial resolution. Recently, iMerenda et al.l (l2006h observed 
He I 10830 A in a polar crown prominence above the limb, 
and found evidence for fields of about 30 G that are oriented 
only 25° from the vertical direction. 
How is the plasma in hedgerow prominences sup- 



ers have suggested that the magnetic field in hedgerow 
prominences is vertical along the observed threads , and 
that the plasma is supported by MHD waves (Jenseiidll 
119861: iPecseli & Engvoldl 120001) . However, relatively high 
frequencies and wave ampUtudes are required, and it is 
unclear why such waves would not lead to strong heating of 
the prominence plasma. Dahlbur g. Antiochos & Klimchukl 
(1998) and Antiochos et al. (1999) showed that the promi- 
nence plasma can be supported by the pressure of a coronal 
plasma lower down along an inclined field line; however, 
this mechan ism only works for nearly horizontal field lines 
(also see Mackav & Galsgaard 2001; Karpenetal. 2005t 
Karpe n. Antiochos & Klimchuk 2006; Karpen & Antiochog 
I2OO8I) . Furthermore, hedgerow prominences are located 
in coronal cavities where the plasma pressure is very low. 
Therefore, it seems unlikely that the vertical threads seen in 
hedgerow prominences can be supported by coronal plasma 
pressure on nearly vertical field lines. 

In this paper we propose that hedgerow prominences are 
embedded in magnetic fields with a complex "tangled" struc- 
ture. Such tangled fields have many dips in the field lines 
where the weight of the prominence plasma can be coun- 
teracted by upward magnetic forces. Our purpose is to 
demonstrate that such tangled fields provide a viable mech- 
anism for prominence support in hedgerow prominences. 
ICasini. Manso Sainz & Low! (|2009|) recently invoked tangled 
fields in the interpretation of spectropolarimetric observations 
of an active region filament. While such filaments are quite 
different from the hedgerow prominences considered here, 
this work shows that tangled fields have important effects on 
the measurement of prominence magnetic fields. Such effects 
will not be considered in this paper 

The paper is organized as follows. In Section 2 we propose 
that hedgerow prominences are supported by tangled mag- 
netic fields, and we discuss how such fields may be formed. 
In Section 3 we present a simple model for the dynamics of 
plasma along the tangled field lines, and we show that weak 
shock waves naturally occur in such plasmas. In Section 4 
we develop a magnetostatic model of tangled fields based on 
the linear force-free field approximation, and in Section 5 we 
study the response of such fields to gravitational forces. In 
Section 6 we simulate the distribution of plasma in a cylindri- 
cal prominence thread. In Section 7 we discuss the formation 
of vertical threads by Rayleigh-Taylor instability. The results 
of the investigation are summarized and discussed in Section 



2. TANGLED FIELDS IN PROMINENCES 

The spectro-polarimetric observations of prominences de- 
scribed in Section [T| are consistent with the idea that quies- 
cent prominences are embedded in coronal flux ropes that 
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Fig. 2. — Magnetic support of solar prominences: (a) by a large-scale 
coronal flux rope, (b) by a tangled magnetic field in a current sheet. 

lie horizontally above the polarity inversion line (PIL). Fig- 
ure lit shows a vertical cross-section through such a flux 
rope. The magnetic field also has a component into the im- 
age plane, so the field lines are helices, and the plasma is 
assumed to be located at the dips of these helical windings. 
A dip is defined as a point where the magnetic field is lo- 
cally horizontal and curved upward. As indicated in the fig- 
ure, the magnetic field may be deformed by the weight of the 
prominence plasma, creating V-shaped dips. The magnetic 
field near such dips is well d escribed by the Kippenhahn- 
Schliiter model (Ki ppenhahn & Schliiter 1957), and several 
authors have developed local magnetostatic models o f the 
fine s t ructures observed in q uiescent promine nces (e.g.,'Lq^ 
Il982t iPetrie & Lowl 120051; iHeinzel. Anz"er & Gunai" 2005^ 
However, recent observations of "dark plumes" ( Berg er et all 
I2OO8) and rotational motions (Chae et al. 2008) within promi- 
nences remind us again that prominences have complex inter- 
nal motions, and it is not clear how such motions can be ex- 
plained in terms of a single large flux rope. Perhaps the mag- 
netic structure of hedgerow prominences is more complicated 
than that predicted by the flux rope model (Figure |2^). 

In this paper we propose an alternative model, 
which is illustrated in Figure Following 
iKuperus & Tandberg-Hansseiil (Il967h . we suggest that 
hedgerow prominences are formed in current sheets that 
overlie certain sections of the PIL on the quiet Sun. Unlike 
those previous authors we suggest that the current sheet 
extends only to limited height 100 Mm), and may extend 
only a Umited distance along the PIL. Furthermore, we 
propose that tangled magnetic fields are present within these 
current sheets. A tangled field is defined as a magnetic 
structure in which the field lines are woven into an intricate 
fabric, and individual field lines follow nearly random paths. 
We suggest that the field is tangled on a spatial scale of 
0.1-1 Mm, comparable to the pressure scale height Hp of the 
prominence plasma {Hp w 0.2 Mm). The prominence plasma 
is assumed to be located at the many dips of the tangled 
field lines. The tangled field is confined horizontally by the 
vertical fields on either side of the sheet, and vertically by the 
weight of the prominence plasma. 

A key feature of a tangled field is that the plasma and field 
are in magnetostatic equilibrium, i.e., the Lorentz force is 
balanced by the gas pressure gradients and gravity. There- 
fore, a tangled field is quite different from "turbulent" mag- 
netic fields in which large-amplitude Alfven waves are present 
(e.g., the solar wind). In a tangled field the magnetic pertur- 
bations do not propagate along the field lines. In this paper 
we examine the basic properties of tangled fields, and we in- 
vestigate their ability to support the prominence plasma. 



We suggest that the tangled field may be formed as a re- 
sult of magnetic reconnection, not the twisting or stress- 
ing of field lines. Quiescent prominences are located 
above polarity inversion "lines" that are often more like 
wide bands of mixed polarity separating regions with dom- 
inantly positive and negative polarity. In these mixed- 
polarity zones, magnetic flux elements move about ran- 
doml y and opposite p olarity elements may cancel each other 
(e.g.. lLivi. Wang & M artin 1985). New magnetic bipoles fre- 
quently emerge from below the photosphere. These pro- 
cesses causes the "recycling" of the photospheric flux about 
once every 2 to 20 hours (Hagenaa r. Schrijver & T itle 2003t 
iHagenaar. DeRosa & Schriiver 20081). a nd the coronal flux is 
recycled even faster (Close et al.l l2005l) . It is likely that the 
interactions between these flux elements produce a complex, 
non-potential magnetic field in the low corona. Within this en- 
vironment magnetic reconnection is likely to occur frequently 
at many different sites in the corona above the inversion zone. 
Each reconnection event may produce a bundle of twisted 
field, and the twisted fields from different events may collect 
into larger conglomerates to form a tangled field. The tangled 
field may rise to larger heights (as a result of its natural buoy- 
ancy), and may collect into a thick sheet that is sandwiched 
between smoother fields, as illustrated in Figure |2j5. The ob- 
served prominence consists of plasma that is trapped within 
this sheet. New tangled field is continually injected into the 
sheet from below, producing vertical motions withi n the sheet. 
We suggest that the "dark plumes" observed by .Berger et al.l 
(I2OO8) may be a manifestation of such vertical motions of the 
tangled field. 

3. FLOWS ALONG THE TANGLED FIELD 

The spatial distribution of plasma within the prominence is 
determined in part by the dynamics of plasma along the tan- 
gled field lines. Figure |3] shows the contorted (but generally 
downward) path of an individual field line in the tangled field. 
Note that there are several "dips" where the field line is hori- 
zontal and curved upward, and "peaks" where the field is hor- 
izontal and curved downward. Tracing the field line upward 
from a dip, one always reaches a peak where the field line 
again turns downward. Therefore, the question arises whether 
the plasma collected in the dips would remain in these dips or 
be siphoned out of the dips via the peaks of the field lines. 

To answer this question, we consider a simple model for the 
motion of the prominence plasma along the magnetic field. 
For simplicity we assume that the flow takes place in a thin 
tube surrounding the selected field line (i.e., the divergence of 
neighboring field lines is neglected), and the cross-sectional 
area of this tube is taken to be constant. We assume a steady 
flow is established along the tube. Let v{s) and p(s) be the 
plasma velocity and density as functions of position s along 
the tube, then conservation of mass requires pv = constant. 
The equation of motion of the plasma is 

dv dp dz 

pv—=--j — Pg-r, (1) 
ds ds 



ds 



where p{s) is the plasma pressure, z{s) is the height above the 
photosphere, and g is the acceleration of gravity. The equation 
of state is written in the form p = Kp^ , where 7 and K are con- 
stants (we use 7 < 5/3 to describe non-adiabatic processes). 
Eliminating p{s) and pis) from equation ([T]l, we obtain the 
following equation for the parallel flow velocity: 

(2) 



dv 
ds' 



^dz 
' ds^ 
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Fig. 3. — Flows along a highly distorted field line in the tangled magnetic 
field of a solar prominence. The field line is indicated by the solid cun'e, 
and the black dots indicate "dips" in the field line where cool prominence 
plasma can collect. The arrows show the direction of subsonic (blue) and 
supersonic {red) flows. Sonic points are located at peaks in the field lines 
{black vertical bars), and shock waves occur where supersonic flows slow 
down before reaching a dip {magenta bars). 

where c{s) is the sound speed (c^ = ^pj p). The above equa- 
tion has a critical point where the flow velocity equals the 
sound speed (v = c). Therefore, a transition from subsonic 
to supersonic flow can occur only at points where the RHS 
of this equation vanishes, dz/ds = 0. These sonic points are 
located at the peaks of the field lines where matter can be si- 
phoned out of one dip and deposited into another dip at lower 
height. The resulting flow pattern is indicated in Figure [3] 
As the supersonic flow approaches the next dip, it must slow 
down to subsonic speeds, which can only happen in a shock. 
Therefore, the tube has a series of subsonic and supersonic 
flows separated by shocks and sonic points. The role of these 
shocks is to dissipate the gravitational energy that is released 
by the falling matter 

The position and strength of the shocks can be computed if 
the height z{s) of the flow tube is known. Neighboring peaks 
are generally not at the same height. Therefore, each sec- 
tion between neighboring peaks is approximated as a large- 
amplitude sinusoidal perturbation superposed on a generally 
downward path: 



z(.s-) « Acos (27r— -00 



-Cs, 



(3) 



where s is the position along the flow tube, A is the distance 
between neighboring peaks (as measured along the flow tube), 
A is the amplitude of the perturbation in height, 00 is a phase 
angle, and C is the background slope. The phase angle is 
chosen such that the peaks in the flow tube (where dz/ds = 0) 
are located at 5 = and s = A, then the slope is given by 
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Fig. 4. — Model for plasma flow along a single field line in a tangled field, 
(a) The solid curve shows the height z(s) as function of position s along the 
field line [in units of the distance A between neighboring peaks, see equation 
(3)]. The dash ed curve shows the overall downward trend of the flow tube, (b) 
The solid cw\e shows the parallel flow velocity v(s) (in units of co) for 7 = 1.5 
and HpQ = 0.15A. The dashed curves show the subsonic and supersonic 
solutions of equation (6). The sonic points are located at the peaks in the field 
line (s = and .v = A), and the shock is located at s = 0.325 A. (c) Similar plot 
for the sound speed c{s) (in units of cq). (d) The solid curve shows the plasma 
density p(s) along the tube (in arbitrary units). The hydrostatic equilibrium 
density is shown by the dashed curve. 

The sonic points will then be located at i = and i = A. Figure 
shows the height z{s) for A = 0.15A and 00 = 0.2 rad, so 
that C = 0.187. Let pa, po and co be the pressure, density and 
sound speed at the sonic points, then K = j 
sound speed can be written as 



'cqPo ^, and the 



Inserting this expression into equation (|2|, we obtain 
du 1 dz 



(u — u 



ds jHp fl ds ' 



(5) 



(6) 



where u{s) = v(s)/co, and Hp.o = Po/ipog) is the pressure 
scale height at the sonic points. Equation (|6]l can be integrated 
as follows: 



I r 2/ \ n , 1 r I-Y/ \ il z(s) — z(So) 

Mm («)-!] + r m' ^{s)-1\ = — 

7-1 jHp 



(7) 



where sq is the position of a sonic point. Equation ^ can be 
solved for u{s) by Newton-Raphson iteration. When 7=1, 
there is an analy tic solution for u(s) in terms of the Lambert 
W function (see ICranmeiil2004h : in the present paper we as- 
sume 7 = 1.5. The supersonic solution vi(.s') is computed with 
io = 0, and the subsonic solution V2{s) is computed with sq = A. 
The dashed curves in Figure|4}5 show vi(.s) and V2(.s-), and Fig- 
ure|4]: shows the corresponding sound speeds ci(s) and C2(s). 
Here we assumed a scale height HpX) = 0.15 A, equal to the 
amplitude of the field-line distortions. The Mach number of 
the flow is given by 



M(s) = v(s)/c(s)=[u{s)] 



(7+I)/2 



(8) 



The shock is located at the point Ssi, where the Mach number 
Ml before the shock and the Mach number M2 after the shock 
satisfy the following relationship: 



C = 27r— sin 00- 



(4) 



2_ 2 + (7- DM? 
27M2-(7-1)' 



(9) 
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which follow s from the Rankine-Hugi onot conditions for par- 
allel shocks (iLandau & Lifshit3ll959h . Therefore, the actual 
flow velocity v(s) between the two sonic points is given by the 
full curve in Figure |4j3, and the sound speed c{s) is given by 
the full curve in Figure |4j;. 

The plasma density p(s) along the flow tube is determined 
by mass conservation (pv = constant), and is plotted in Figure 
|4}l. Note that there is a strong peak in the density at the dip 
in the field line, i^,-^ « 0.56A. The dashed curve in Figure|4jl 
shows the density profile that would exist if the plasma were 
in hydrostatic equilibrium (HE): 



Phe(s) = Pdip^W 



z(s)-z(sdip) 



(10) 



where pdip and Hp are the density and pressure scale height at 
the dip. Hp = 0.254A. The deviations from hydrostatic equi- 
librium are significant only in those regions where the flow 
velocity is comparable to the sound speed. We define an aver- 
age flow velocity v by 



A 



(11) 



For the case shown in Figured we find v = 0.6co, so the av- 
erage flow speed is less than the sound speed. Therefore, the 
contorted shape of the flow tube significantly reduces the flow 
velocity compared to the supersonic free fall that would occur 
in a straight vertical tube. 

The cooler parts of the prominence are thought to have a 
temperature T ^ 10^ K. Assuming a hydrogen ionization frac- 
tion of 10%, a helium abundance of 10 % and 7 = 1.5, the 
sound speed cq w 10 km s"', and we predict an average flow 
velocity v « 0.6co ~ 6 km s"' . The vertical component of this 
velocity is v. -Cv w -1.1 km s~', less than the observed 
vertical velocities in prominence threads (5-10 km s"'). Note 
that the predicted velocity is relative to the pattern of the tan- 
gled field, therefore, if the tangled field expands in the ver- 
tical direction it will push the prominence plasma upward. 
We speculate that t he observed upwar d motions in hedgerow 
prominences (e.g.. iBerger et al.ll2008h are due to such large- 
scale changes in the tangled field. 

4. LINEAR FORCE-FREE FIELD MODELS 

In this section simple models for tangled fields are devel- 
oped. A volume V in the corona is considered, and the plasma 
inside this volume is assumed to be in magnetostatic equilib- 
rium, -Vp+pg+F = 0, where p is the plasma pressure, p is the 
density, g is the acceleration of gravity, and F is the Lorentz 
force. All quantities are functions of position r within the 
volume. The Lorentz force is given by 



F=ijxB=— (VxB)xB, 

C A-TT 



(12) 



where j is the electric current density and B is the magnetic 
field. Using tensor notation, equation ( fT2b can also be written 
as Fi = dTij/dxj, where Tij is the magnetic stress tensor, a 
special case of Maxwell's stress tensor ( Jackson. 1999) : 



?1 



-5ii + 



47r 



(13) 



The first term describes magnetic pressure, and the second 
term describes magnetic tension. In a tangled field both pres- 
sure and tension forces are important. 



If gravity and plasma pressure gradients are neglected, then 
F 0, so the magnetic field B(r) must satisfy the force-free 
condition: 

VxB = aB, (14) 

where a(r) may in general be a function of position. In the 
special case that a is constant throughout the volume, equa- 
tion ( fT4b becomes a linear equation for B(r) , and the solu- 
tions are called linear force free fields (LFFF). IWoltieij ( Il958h 
has shown that in a closed magnetic system with a prescribed 
magnetic helicity (H = J A- BdV, where A is the vector po- 
tential) the lowest-energy state is a LFFF. Therefore, in this 
paper only LFFFs are considered, and a is treated as a free 
parameter. We find that LFFFs can be tangled. The typical 
length scale of the tangled field is given by the inverse of the 
a parameter, £ = |q;|~'. In the following we first consider the 
case that £ is small compared to the domain size L in any di- 
rection, and then consider the boundary effects. Section |43] 
describes tangled fields in a cylindrical domain. 

4.1. Tangled Field in a Large Volume 

In the absense of boundary conditions, the solution of equa- 
tion (O can be written as a superposition of planar modes: 

N 

B(r) = ^B„ [ei,„cos(k„-r+/3„)-e2,„sin(k„-r+/3„)] , (15) 

n=l 

where is the number of modes, k„ = acs.,, is the wave vector 
(n = 1 , • • • , A^), B„ is the mode amplitude, /3„ is a phase angle, 
and [ei,„,e2,«,e3,„] are unit vectors that are mutually orthogo- 
nal and form a right-handed basis system: 

(16) 
(17) 
(18) 



Ci.B = cos6'„(cos(/)„ f + siiKpn z)-sm9„ x, 

e2.«=-sin0„ y+cos0„ z, 

e3.„ = sin0„(cos0„ y + sin0„ z) + cos6'„x. 



Here 6„ and (j)„ are the direction angles of the wave vector 
relative to the Cartesian reference frame. 

Figure |5] shows an example of a field with = 100 modes, 
an isotropic distribution of direction angles (6*,,, </)„), and ran- 
domly selected phase angles /3„. The starting points of the 
field lines are randomly selected from the central part of the 
box, and the field lines are traced until they reach the box 
walls. Note that individual field lines follow random paths, 
and that different field lines are tangled together 

Now consider an ensemble £^jv of fields with different 
(randomly distributed) phase angles /?„, but with a fixed num- 
ber of modes A^, and with fixed mode amplitudes B„ and di- 
rection angles (6*,,, </)„). The phase angles are assumed to be 
uniformly distributed in the range [0, In], and angles from dif- 
ferent modes n and n' are assumed to be uncorrected. Then 
the ensemble average of the magnetic field vanishes. 



<B>^=0, 



(19) 



where < ••• >^ denote the average over phase angles f3„. 
Also, the average of the tensor BB is given by 

N 

<BB>;3= 5^B^(ei,„ei,„ + e2,„e2,„) 



«=i 



= 5^B^(l-e3.„e3,„) 



(20) 
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Fig. 5. — Tangled field obtained by superposition of 100 randomly selected 
modes of the linear force-free field in a large volume. 

where I is the unit tensor (I = xx+yy+zz). Expression ( l20b is 
independent of position r, so the magnetic field is statistically 
homogeneous, but is not necessarily isotropic. 

We study the statistical properties of field lines in models 
with different values of the mode number A^. For each 
we construct a series of models (m = 1, • • • ,M) with differ- 
ent phase angles, but with constant values of the mode ampli- 
tudes B„ and direction angles (d„,(j)„). The mode amplitudes 
are taken to be the same for all modes {B„ = 1). For each real- 
ization m of the phase angles (3„ we trace out the field line that 
starts at the origin {x = y = z = Q), and we measure the square 
of the radial distance rl,is) as function of position s along the 
field line: 

ri(s)=xl(s)+yl(s) + zi(s). (21) 

We then average this quantity overM = 100 realizations of the 
phase angles to obtain the mean square distance n(s). ForN = 
3 both mutually orthogonal directions {d„,(f)„) and randomly 
chosen directions are considered. In both cases^ we find that 
the field lines follow long helical paths, and r^{s) increases 
quadratically with s. Therefore, for = 3 the field lines do not 
behave randomly. For > 4 only randomly chosen directions 
are considered. For N = 4 some of the field lines are long 
helices, while others have more random paths, and for A^ = 
5 all field lines seem random, however, in both cases r^(s) 
is not well fit by a power law. True random walk behavior 
of the field lines, as indicated by a linear dependence of p- 
on distance s, is found only when the number of modes is 
increased to N > 10. In the limit of large A^, r^(s) sa lOi/laj. 

We now consider a larger ensemble £n in which not only the 
phase angles /3„ but also the mode amplitudes B„ and direction 
angles i9n,(j)„) are allowed to vary. From now on < • • • > will 
denote an average over this larger ensemble. The direction an- 
gles are assumed to have an isotropic distribution, i.e., the an- 
gle (/)„ is uniformly distributed in the range [0,27r], and cos 6',, 
is uniform in the range [-1 ,+1], so that < cos^ 6„ >= |. Fur- 
thermore, the mode amplitudes B„ are assumed to be uncorre- 



lated with the direction angles, and < B^ >= B^/N, where Bo 
is a constant. Then < B >= 0, so the mean magnetic field van- 
ishes. Further averaging of equation (l20l l shows that < BB > 
is an isotropic tensor: 

N 

< BB >= i ^ < > (i_ < e3 „e3„ >) = \bIi. (22) 

«=i 

It follows that < B^ >= Bq, so Bo equals the r.m.s. value of 
the total field strength. The ensemble average of the magnetic 
stress tensor, equation ( fT3] ), is given by 



which is also isotropic. Note that the diagonal components 
of < Tij > are negative, so the effects of magnetic pressure 
dominate over the effects of magnetic tension. Therefore, the 
isotropic tangled field has a positive magnetic pressure, p, = 
B\/{2Ati). The average energy density is E, = bIJ{%-k). The 
relati onship E, = 3p, is similar to that for a relativistic gas 
(e.g.. lWeinber"g[T972h . 

4.2. Boundary Effects 

The tangled field must be confined within a certain vol- 
ume (e.g., a current sheet, see Figure Wp), and the confine- 
ment must be effective for a period much longer than the 
Alfven travel time across the volume. What are the condi- 
tions for such confinement? To answer this question we must 
consider the boundary region between a tangled field and a 
smooth field. The tangled field is assumed to be character- 
ized by a high value of |a|, and the smooth field presumably 
has a much lower value of To last a long time, the mag- 
netic field near the boundary must be approximately in equi- 
librium (non-linear force-free field). The force-free condition 
(O implies that a is constant along field lines, so there can- 
not be many field lines that pass from the smooth region to 
the tangled region. Therefore, one important condition for the 
survival of the tangled field is that the two regions are nearly 
disconnected from each other magnetically. Another require- 
ment is that the two regions are approximately in pressure bal- 
ance. 

To show that these conditions can be satisfied, we now con- 
sider a simple model for the boundary region. The interface 
between the tangled and smooth fields is approximated by a 
plane surface, here taken to be the plane x = Om Cartesian co- 
ordinates. The above-mentioned condition on the lack of con- 
nectivity between the smooth and the tangled fields requires 
Bjc(Q,y,z) = at X = 0. The tangled field in jc > is assumed 
to be a LFFF with a specified value of a. The solution of the 
LFFF equation is again written as a superposition of planar 
modes. However, in the present case the modes are grouped 
into pairs with closely related wave vectors k„ and k^, and 
with the same amplitude B„ and phase /?„: 

N/2 

B(r) = ^B„ [ei.„cos(k„ •r + ^„)-e2,„sin(k„ •r + /?„) 
«=i 

e'i „cos(k;-r + ;3„) + 62,,, sin(k;,-r + /?„)] . (24) 

Here A^ is the total number of modes, ei ,„ and ea.n are defined 
in equations (fTSI l and ([TtT i. k^ = ae'j^ „ is the modified wave 
vector, and the unit vectors e'j „ and §3 „ are defined by 

e'l „=-cos0„(cos(/)„ y + sin(/)„ z)-sin0„ x, (25) 

e3„= sin0„(cos0„ y + sin0„ z)-cos0„ X. (26) 
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Note that 63 „ differs from e3.„ only in the sign of 
the jc-component, whereas e'j „ has the sign of the y- 
and z-components reversed. Therefore, the unit vectors 
[e'j „,e2,„,e3 „] again form a right-handed basis system. The 
magnetic field at the boundary (x = 0) is given by 

N/2 

B„cos(?„(cos(/!>„ y + sm(j)„ z) 



X cos [Q;sin6'„(ycos0„ + zsin(/)„) + /3„] , 



(27) 



which satisfies BJO,y,z) = 0. Therefore, it is possible to con- 
struct a tangled field that is disconnected from its surround- 
ings. 

We now consider the statistical average of the tensor BB at 
x = 0. Averaging over phase angles /3„, we obtain 

N/2 

<BB >f3 = 2 bI cos^ On (cos (p„ y + sm(p„ z) 



X (cos0„ y+sin0„ z) atx = 0, 



(28) 



and further averaging over mode amplitudes and direction an- 
gles yields 

<BB >= iB^(yy + zz) atjc = 0. (29) 
6 

Here we assume an isotropic distribution of direction angles, 
and we use < >= Bf^/N, where Bq is the r.m.s. field strength 
in the interior of the tangled field (see Section im i. Note that 
at the boundary < B^(0,y,z) >= bI/3, while in the interior 
< B^ >= Bq, so the rm.s. field strength at the boundary is 
reduced by a factor 1 / a/3 compared to that in the interior. 
The magnetic pressure at x = is given by 

<B\0,y,z)> Bl 

5 = ^^^^ 

OTT 247r 

where pi is the average pressure in the interior of the tangled 
field [see equation (|23]|1. Equation ( l30l l shows that it is pos- 
sible to maintain pressure balance between the tangled field 
and its surroundings. 

4.3. Tangled Field in a Cylinder 

Here an infinitely long cylinder with radius R is consid- 
ered. We adopt a cylindrical coordinate system (r,(j),z), and 
we assume that the radial component of magnetic field van- 
ishes at the cylinder wall, Br{R,4),z) = 0. In the Appendix 
we analyze the eigenmodes of the LFFF equation in the do- 
main r <R subject to the above boundary condition. We find 
that this eigenvalue problem has a discrete set of modes, and 
the number of modes depends on the dimensionless param- 
eter a = \a\R. Figure |6] shows the resulting magnetic fields 
for fl = 3.0, 4.5 and 6.0. In the first case only the axisymmet- 
ric (Lundquist) mode is present, so the field lines are helical. 
Assuming the cylinder axis is vertical, there are no dips in the 
field lines. If cool plasma were to be injected into such a struc- 
ture, it would spiral down along the field lines and quickly 
reach supersonic speeds. In contrast, for a = 4.5 and a = 6.0 
there are multiple modes of the LFFF, and the random super- 
position of these modes creates a tangled field with many dips 
where prominence plasma can be supported. The field-line 
dips (i.e., sites where B, = and B • VB^ > 0) are indicated by 
magenta dots in the middle and right panels of Figure |6] We 
will return to this model in Section|2] 






(a) 



(b) 



(c) 



Fig. 6. — Tangled magnetic fields obtained by superposition of modes of 
the linear force-free field inside a cylinder, which represents a vertical thi'ead 
within a hedgerow prominence. The radial component of the field B, vanishes 
at the cylinder wall. The three panels show models with different values of 
a = \a\R, where a is the torsion parameter and R is the cylinder radius: (a) 
a = 3.0; (b) a = 4.5 ; (c) a = 6.0. In case (a) the field has only a single mode, 
but is cases (b) and (c) there are multiple modes {N = A and N = 6), some 
of which are non-axisymmetric, resulting in tangled field lines. Each panel 
shows three field lines {red, green and blue curves). In panels (b) and (c) the 
magenta dots show dips in the field lines. 

5. DEVIATIONS FROM THE FORCE-FREE CONDITION 

The above models for a tangled field are purely force-free 
and do not have any magnetic forces to support the promi- 
nence plasma against gravity. To include such effects, we 
now consider the "elastic" properties of the tangled field, i.e., 
its response to external forces. Specifically, the weight of the 
prominence causes the tangled field to be compressed in the 
vertical direction, resulting in a radially outward force on the 
plasma. Also, shearing motions may occur within the tan- 
gled field as dense plasma m oves downward an d less dense 
"plumes" move upward (e.g., Berger et"ani2008h . This results 
in shear deformation of the tangled field and associated mag- 
netic stresses that counteract the plasma flows. In the follow- 
ing both of these effects are considered in some detail. 

5.1. Compressional Effect 

We first consider the effects of gravity on a layer of tangled 
magnetic field. The magnetostatic equation (-Vp+pg+F = 0) 
cannot be solved analytically for a tangled field, so we make 
the following approximation: 



B'(r) = V X 



-B(r)e 

a 



-z/Hb 



B(r)- 



1 



aHB 



z X B(r) 



-z/Hb 



(31) 



Here B(r) is the LFFF given by equation ( fTSl ), and Hb is 
the magnetic scale height of the modified field (we assume 
\a\HB > 2). Note that V • B' = as required. This modified 
field B'(r) is no longer force-free, but has the following sta- 
tistical properties: 

< B' >= 0, (32) 

< B'B' >= \bI [{l + e\xiL+yy) + ii\ e-^'^"" , (33) 
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where e 



\a\HBr' < ' 



2- 



Therefore, the magnitude of the 



modified field drops off exponentially with height z- Let T/j 
be the magnetic stress tensor of the modified field. Taking its 
statistical average, we find for the nonzero components of the 
stress tensor: 



<r;,>=<r; >=S- ^'^'""^ 

247r 

24tt 



(34) 



(35) 



Note that for e < j the stress tensor is nearly isotropic. The 
net force on the plasma is given by F/ = dT/j/Oxj, and the 
average force follows from equations ( l34l i and ( l35T i: 



<f^>=<f;. >=o, 



<f!>-- 



Bi 



UttHb 



(l+2e^) e 



(36) 
(37) 



Note that the average force acts in the positive z direction, i.e., 
the magnetic force counteracts the force of gravity. In effect, 
the plasma is being supported by the magnetic pressure of the 
tangled field. The tangled field acts like a hot gas that has a 
significant pressure but no mass. The average density of the 
plasma that can be supported by the tangled field is given by 



Pavgiz) ■■ 



<f: > 



( l + 2e^) e 



-Iz/Hb 



(38) 



where g is the acceleration of gravity. 

The horizontal components of Lorentz force, F'^ and F,', do 
not vanish for any particular realization of the tangled iSeld, 
and cannot be written as the gradients of a scalar pressure p. 
The reason is that expression (ISTT l is not an exact solution of 
the magnetostatic equilibrium equation. However, equation 
(|36] ) shows that the horizontal forces vanish when averaged 
over the fluctuations of the isotropic tangled field. Therefore, 
expression (ISTT i is thought to give a good approximation for 
the effects of gravity on the tangled field. 

We now apply the above model to the vertical threads ob- 
served in hedgerow prominences. To explain the observed 
heights of such prominences, we requke that the magnetic 
scale height Hb is at least 100 Mm. The size ^ of the mag- 
netic tangles is assumed to be in the range 0.1-1 Mm, so 
e = i/HB < 1. For Bq = 10 G and Hb = 100 Mm, we find 
Pavg ~ 10~'^ g cm~^, which corresponds to an average (total) 
hydrogen density nn.avg ~ 5 x lO' cm"^. This is only about 
0.05 times the density po ~ 2 x 10"'^ g cm~^ or «^ « 10'^ 
cm~^ typically observed in h edgerow prominences (Engvold 
1 1 97611 1 98d: iHirayamal 1 9861) . This comparison shows that the 
pressure of the tangled field inside a prominence thread is not 
sufficient to support the weight of the prominence plasma. To 
support the plasma with tangled fields, we need to take into 
account the magnetic coupling between the vertical thread 
and its surroudings. Such coupling is neglected in the above 
plane-parallel model. 

5.2. Shear Stress Effect 

We now assume that the tangled field pervades not only 
the observed vertical threads but also their local surroundings. 
The density in the surroundings is less than that in the threads, 
so the force of gravity is also much lower. This difference 
in gravitational forces leads to vertical motions (downflows 
in the dense threads, upflows in the tenuous surroundings) 



that create magnetic stresses in the tangled field. The mag- 
netic coupling between the prominence and its surroundings 
causes the weight of the dense prominence to be distributed 
over a wider area. In effect, the prominence plasma is being 
supported by the radial gradient of the magnetic pressure of 
the tangled field over this larger area. In the following we 
estimate the magnetic stresses and vertical displacements re- 
sulting from these forces. 

The tangled field is modeled either as a vertical slab with 
half-width R, or as a vertical cylinder with radius R. The 
prominence is located at the center of this slab or cylinder, 
and has a half- width or radius ro < R. Then the average den- 
sity in the tangled field region is given by 



Pavg ~ PO 



(39) 



where po is the density of the prominence, and we neglect the 
mass of the surroundings. The exponent n = 1 for the slab 
model or n = 2 for the cylindrical model. As discussed in 
Section 15.11 observations of hedgerow prominences indicate 
po ~ 2 X lO"'-' g cm"^ (e.g., Engvold 1976), and to explain the 
observed heights of such prominences with Bo= 10 G, we re- 
quire Pavg/ Po < 0.05. According to equation ( [39] l. this implies 
R/rQ>2Q for the slab model, or R/ro > 4.5 for the cylindri- 
cal mod el. The observe d threads have widths down to about 
500 km dEngvoldll 19761) . which corresponds to ro w 250 km. 
Therefore, the magnetic coupling by the tangled field must 
extend to a surrounding distance of at least 5 Mm for the slab 
model, or 1.1 Mm for the cylindrical model. More generally, 
equations ( |38] ) and ( [39] ) yield the following expression for the 
magnetic scale height of tangled field: 



Hb 



R 



(40) 



Therefore, the maximum height of the prominence depends 
strongly on the magnetic field strength Bq. 

According to the present model, magnetic stress builds up 
in the tangled field as a result of the difference in gravitational 
forces between the thread and its surroundings. Can the field 
support such shear stress? To answer this question we exam- 
ine the effect of vertical displacements on the tangled field. 
For simplicity we neglect the mean vertical force given by 
equation (|37] |. and we focus on relative displacements. Let 
r' be the new position of a fluid parcel originally located at 
position r. In the limit of a perfectly conducting plasma, the 
deformed field B' at the new position r' is given by 



B' = 



B^d^^ 



(41) 



where B(r) is the original field, and J is the Jacobian of the 
transformation (e.g.. Priest 1982). We assume 



-■X, y' =y, z' = z+h(x,y), 



which yields 

B'^ = B^, B[.=By, 



, dh dh 

Bl = B: + B,— +By — , 

ox ' oy 



(42) 



(43) 



where h(x,y) is the vertical displacement. The original field 
B(r) is assumed to be a realization of the isotropic tangled 
field given by equation (fTSI) . Using equation (|22] |. we obtain 

< B' >= 0, 



< B'B' >= 



1 r2 



9/!.: 

XX + yy + /zz + — xz + — yz 

ox oy 



(44) 
(45) 
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where 



/(^,3')=1 + 



(46) 



This yields the following expressions for the off-diagonal 
components of the stress tensor: 



Bl dh , Bl dh 
< j' >= — 2 < j' >= — y 



(47) 



The Lorentz force is given by Fj' = dTlj/dxj, and since < T/, > 
is independent of z, the average vertical force is given by 



Bl fd^h d^h\ , , , 



(48) 



where /S.p{x^y) is the density perturbation (Ap = p-pavg)- Us- 
ing this equation, we can determine the vertical displacement 
h{x,y) for a given density variation Ap{x,y). In the follow- 
ing subsections we solve the above equation for the slab and 
cyhnder models. 

5.2.1. Slab Model 

We first consider a slab with infinite extent in the +y, -y and 
+z directions. The coordinate x perpendicular to the slab is in 
the range -R < x < R, where R is the half-width of the sheet 
in which the tangled field is embedded (see Figure|2]3). Then 
the density perturbation is given by 



Ap(x)=(+^''[V('-»/^)] 
-Paro R 



for \x\ < ro, 
otherwise. 



(49) 



Inserting this expression into equation ( |48] l and solving for the 
vertical displacement, we obtain 

C[l-(ro/R)](x^-rl) 
C(ro/R)[(R-rof-iR- 



h(x)-- 



for |x| < ro, 
jjcj)^] otherwise, 



(50) 



where C = 6TTgpoB'[^. Here we applied no-stress boundary 
conditions (dh /dx=Q)atx= ±R. Note that h{x) and its deriva- 
tive are continuous at the edges of the prominence {x = ±ro). 
The relative displacement across the tangled field is given by 



Ah = h(R)-h(0)-- 



67rgpo 



ro(R-rQ). 



(51) 



Figure shows the function h(x) for ro = 0.5 Mm, R = 10 
Mm, Bo = 10 G and po = 2x 10"'^ g cm-^ so that Ah = 0.491 
Mm. Note that a relatively small deformation of the tangled 
field {Ah ^ R) is sufficient to redistribute the gravitational 
forces over the full width of the tangled field. However, Ah is 
larger than the pressure scale height of the prominence plasma 
(Hp « 0.2 Mm). Therefore, the deformation of the tangled 
field in the neighborhood of the prominence may have a sig- 
nificant effect on the distribution of the prominence plasma. 
This issue will be further discussed in Section |7] 

For comparison of the tangled field slab model with the flux 
rope model (Figure |2|i, we define the average sag angle 9 of 
the prominence relative to its surroundings: 



tanfy = 



Ah 6Trgporo 



R 



Bi 



(52) 



where Bo is the r.m.s. field strength of the tangled field, and 
we assumed ro <C R- This expression is similar to that derived 
for the Kippenhahn-Schliiter model: 



tan6' = 



Bl ' 



(53) 



(a) 




(b) 




Fig. 7. — Vertical displacements in the tangled magnetic field supporting 
a hedgerow prominence, (a) Slab model with R = 10 Mm, = 0.5 Mm, and 
Bo = 10 G. The vertical displacement h(x) is plotted as function of position ,v 
perpendicular to the vertical slab. The location of the prominence is indicated 
by the dotted lines, (b) Cylindrical model of a prominence thi'ead with R = 5 
Mm. ro = 1 Mm, and Bo = 3 G. In both cases the prominence density po = 
2x gcm-^ 



where B^ is the magnetic field through the mid-plane of the 
prominence (Kippenhahn & Schluter 1957). Equation ( l53T l 
describes the angle of the field lines in the flux rope model 
shown in Figure|2^. Therefore, the flux-rope and tangled field 
models are similar in their ability to explain the magnetic sup- 
port of the prominence plasma, provided the half-width R of 
the tangled field region is similar to the radius R of the flux 
rope. 

5.2.2. Cylindrical Model 

We now consider a cylindrical model for a prominence 
thread with r the distance from the (vertical) thread axis. The 
density perturbation is given by 



Ap(r) = 



r +pa[\-{n/Rf] forr<ro, 
\ -poiro/Rf for ro<r<R. 



(54) 
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The vertical displacement is obtained by solving equation 
(081), which yields 



h(r): 



where 



/i„[(r/ro)2-l] 



for r < ro, 



ho[2RHn(r/ro)-rHrl]/(R^-rl) for r,, < r < R, 



(55) 
(56) 



and we applied no-stress boundary conditions at r = R. Then 
the total displacement across the tangled field is 



Ah = h(R)-h(Q) = 



T — ' 



;ln(^ 

^0 



(57) 



Figure|7J' shows the vertical displacement h(r) for R = 5 Mm, 
ro = 1 Mm, Bo = 3 G and po = 2x 10"'-^ g cm"-', which yields 
ho = 0.55 1 Mm and Ah = 1 .847 Mm. In this case Ah is signif- 
icantly larger than the pressure scale height (Hp ss 0.2 Mm), 
mainly because of the lower field strength compared to the 
case shown in Figure |7^. In Section |7] we consider the effect 
of such deformation on the field-line dips, and on the spatial 
distribution of the prominence plasma. 

The above analyses only provide an rough estimate for the 
density of prominence plasma that can be supported by the 
tangled field. The actual density distribution p(r,f) is likely to 
be much more complex for several reasons. First, plasma will 
tend to collect at the dips of the field lines, so the density will 
vary on the spatial scale £ of the tangled field and on the scale 
of Hj,; this effect will be considered in more detail in Section 
|7] Second, the density will vary with time because there are 
flows along the field lines (see Section O and these flows are 
likely to be non-steady. Also, the magnetic structure is not 
fixed and will continually evolve as dipped field lines are dis- 
torted by the weight of the prominence plasma. To predict 
the actual density will require numerical simulations of the 
interaction of tangled fields with prominence plasma, which 
is beyond the scope of the present paper. 

6. FORMATION OF VERTICAL THREADS BY RAYLEIGH-TAYLOR 
INSTABILITY 

According to the present theory, hedgerow prominences are 
supported by the pressure of a tangled magnetic field, which 
acts like a tenuous gas and is naturally buoyant. It is well 
known that a tenuous medium supporting a dense medium 
is s ubject to Rayleigh-Taylor (RT) instability ( Chandrasekh^ 
1196 lb . Therefore, we suggest that the observed vertical 
threads may be a consequence of RT instability acting on the 
tangled field and the plasma contained within it. As cool 
plasma collects in certain regions of the tangled field, the 
weight of the plasma deforms the surrounding field, which 
causes even more plasma to flow into these regions. 

A detailed analysis of the formation of prominence threads 
by RT instability is complicated by the fact that we presently 
do not understand how a tangled field responds to shear defor- 
mation. In Section I572I we estimated the relative vertical dis- 
placement Ah of the prominence plasma assuming no recon- 
nection occurs during the deformation of the magnetic field 
by gravity forces [see equations ( BTT i and (|57]i1. In this case 
the tangled field behaves as an "elastic" medium with mag- 
netic forces proportional to the displacement. However, it 
is not clear that this approximation is valid. High-resolution 
observations of prominences indicate that the dense threads 
move downward relative to their more tenuous surroundings 



with velocities of the or der of 10-30 km s~' (e.g..l Berger et al.l 
l2008HChae et al.ll2008h . If the threads and their surroundings 
are indeed coupled via tangled fields, these relative motions 
imply that the field is continually being stretched in the verti- 
cal direction. Therefore, the shear stress continually increases 
with time, unless there is internal reconnection that causes the 
shear stress to be reduced. 

We speculate that tangled fields have a tendency to relax to 
the LFFF via internal reconnection. A similar relaxation pro- 
cesses occurs in the reve rsed field pinc h and other laboratory 
plasma physics devices (|Ta¥lor||1974l). Therefore, the long- 
term evolution of prominence threads likely involves small- 
scale reconnection within the tangled field. The tangled field 
may behave more like a "plastic" medium that is irreversibly 
deformed when subjected to shear stress. Such plasticity 
makes it possible to understand how the dense threads can 
move downward relative their the surroundings at a small but 
constant speed. These flows significantly deform the tangled 
field, but the field is nevertheless able to support the plasma 
against gravity. A detailed analysis of reconnection in tangled 
fields and its effect on the prominence plasma is beyond the 
scope of the present paper 

The observed vertical structures likely reflect the non-linear 
development of the RT instability in hedgerow prominences. 
To establish a vertical column of mass resembling a promi- 
nence thread will likely require vertical motions over a sig- 
nificant height range (tens of Mm). Starting from a homoge- 
neous density distribution, it may take several hours for the 
threads to form by RT instability. 

7. MODEL FOR A PROMINENCE THREAD 

We now construct a model for the density distribution in a 
fully formed (vertical) prominence thread supported by a tan- 
gled field. It is assumed that the RT instability has produced 
a vertical thread that is clearly separated from the rest of the 
prominence plasma. Therefore, only a single thread and its 
local surroundings are considered, and the tangled field is as- 
sumed to be contained in a vertical cylinder with radius R = 5 
Mm. As discussed in Section[3] there will in general be mass 
flows along the tangled field lines, but for the purpose of the 
present model we neglect such flows and we assume that the 
plasma is in hydrostatic equilibrium along the field lines. 

To construct the density model, we first compute a partic- 
ular realization of the LFFF with aR = 9 (see Appendix for 
details). To account for the weight of the prominence plasma, 
this field is further deformed as described by equation ( l55b . 
The deformation parameters are ro = 1.25 Mm and ho =1.5 
Mm, which yields Ah = 4.44 Mm, somewhat larger than the 
values used in Figure|7J). As shown in Section lSTl the weight 
of the prominence plasma causes the strength of the tangled 
field to decrease with height [see equation ([3TTi1. but for sim- 
plicity this gradient is neglected here. The pressure scale 
height is assumed to be constant. Hp = 0.2 Mm, which cor- 
responds to a temperature of about 8000 K, typical for Ha 
emitting plasma in prominences. 

We introduce cartesian coordinates (x,y,z) with the z axis 
along the cylinder axis; the x and y coordinates are in the 
range [-R,+R], and the height z is in the range [0, lOR]. The 
density p(x,y,z) in this volume is computed on a grid with 
200 X 200 X 1000 grid points, using the following method. We 
randomly select a large number of points within the cylinder 
and trace out the field lines that pass through these points. For 
each field line we plot the height z{s) as function of position 
s along the field line, and we find the peaks and dips in the 
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Fig. 8. — Height z{s) as function of position i along a field line in the cylin- 
drical tangled field model with a = 9.0. Lengths are in units of the cylinder 
radius R. 

field line. For each dip we find the two neighboring peaks (si 
and S2) and we determine the depth Az of the valley between 
these peaks. We then iteratively remove shallow dips with 
Az < 3Hp by concatenating neighboring sections. Figure |8] 
shows an example for one particular field line; the remaining 
dips with Az > 3Hp are indicated by squares. We then com- 
pute the density p(s) in each interval, assuming hydrostatic 
equilibrium along the field line [see equation (fTOlll. Since we 
are interested only in relative densities, we set pdip = 1, the 
same for all dips on all field lines. Finally, we distribute the 
density p{s) onto the 3D grid by finding the grid points that lie 
closest to the path of the field line. This process is repeated 
for 8000 field lines to obtain the density p(x,y,z) throughout 
the 3D grid. 

Figure|9]shows the resulting density distribution. The three 
panels show different projections obtained by integration in 
the X, y and z directions, respectively [for example. Figure [9^ 
shows J p{x,y,z)dx]. Note that the plasma is concentrated 
in the central part of the cylinder; this is due to the defor- 
mation of the magnetic field described by the vertical dis- 
placement h(r), which changes the distribution of the field- 
line dips compared to that in the LFFF The plasma is highly 
inhomogenous {pmax ~ 50payg), and is distributed in sheets 
corresponding to surfaces of field-line dips. In some regions 
there are multiple sheets along the line of sight (LOS). This is 
consistent with observations of the H I Lyman lines in solar 
prominences, wh ich indicate multiple threads along the LOS 
(lOrrall & Schma hl 1980; Gunaret al. 2007). 

The above model is quasi-static and does not take into ac- 
count the expected mass flows along the field lines (Section 
O, nor the dynamical processes by which the threads are 
formed (Section|6]l. Constructing a more realistic 3D dynam- 
ical model of a prominence thread supported by tangled fields 
is beyond the scope of the present paper However, our main 
conclusion that the density distribution within the thread is 
highly inhomogeneous is likely to be valid also in the dynam- 
ical case. 

8. SUMMARY AND DISCUSSION 

We propose that hedgerow prominences are supported by 
magnetic fields that are "tangled" on a spatial scale of 1 Mm 
or less. A key feature of the model is that the plasma is ap- 
proximately in magnetostatic balance, therefore, the model 
is different from earlier models in whi ch the plasma i s sup- 
ported by MHD-wave p ressure (e.g., iJensenI I1983L Il986t 
iPecseli & Engvoldll2000l) . In the present case the perturba- 
tions of the magnetic field do not propagate along the field 
lines. The tangled field is located within a large-scale current 
sheet standing vertically above the PIL (FigurelJt), and is not 
magnetically connected to the photosphere on either side of 
the PIL. Such tangled fields may be formed by flux emergence 




(q) (b) (c) 

Fig. 9. — Simulated density distribution of a prominence thread. The 
thread is supported by a tangled magnetic field with R = 5 Mm and aR = 9. 
The field is distorted by the weight of the prominence plasma (ro = 1 .25 Mm. 
/iQ = 1 .5 Mm). The plasma is concentrated at the dips of the tangled field lines. 
The density distribution along the field lines assumes hydrostatic equilibrium 
with Hp = 0.2 Mm. The three panels show the column density distributions 
as seen from the x, y and z. directions, respectively. Note that multiple sheets 
of plasma are superposed along the line of sight. 

followed by magnetic reconnection in the low corona. 

In this paper we use a variety of methods to explore the in- 
teractions of prominence plasma with tangled fields. In Sec- 
tion |3] a simple model for the downward flow of plasma along 
the distorted field lines was developed. It was shown that such 
flows naturally develop standing shocks where the gravita- 
tional energy of the plasma is converted into heat; this may be 
important for understanding the heating of prominence plas- 
mas. The average flow velocity is less than the sound speed, 
indicating that the tangled field is able to support the promi- 
nence plasma against gravity. 

In Section|4]linear force-free models of tangled fields were 
constructed. Tangled fields can be described as a superposi- 
tion of planar modes. We studied the statistical properties of 
such fields, and found random-walk behavior of the field lines 
when the number of modes is sufficiently large (A^ > 10). To 
produce tangling of the field lines on a scale of 1 Mm or less, 
as required for our model of hedgerow prominences, we need 
|q!| > 1 Mm"', much larger than the values typically found in 
measurements of photospheric vector fields. Therefore, ac- 
cording to the present model, hedgerow prominences are em- 
bedded in magnetic fields with high magnetic helicity den- 
sity. We also considered the conditions for confinement of a 
tangled field by the surrounding smooth fields, and showed 
that despite the high heUcity density tangled fields can be in 
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pressure balance with their surroundings. 

In Section|5]the "elastic" properties of a tangled field were 
considered, i.e., their linear response to gravitational forces 
assuming ideal MHD. We found that the weight of the promi- 
nence plasma can be supported by the nearly isotropic mag- 
netic pressure of the tangled field. The tangled field pervades 
not only the observed vertical threads, but also their local sur- 
roundings. The magnetic coupling between the threads and 
their surroundings is quite strong: vertical displacements of 
only 0.5-2 Mm are sufficient to counteract the shear stress re- 
sulting from the different gravitational forces. In effect, the 
weight of a dense thread is distributed over an area that is 
larger than the cross-sectional area of the thread. As discussed 
in Section 152] the observed densities in prominence threads 
10" cm"^) can be supported by tangled fields with field 
strengths in the range 3-15 G. 

In Section |6] we proposed that the observed vertical 
structures in hedgerow prominences are a consequence of 
Rayleigh-Taylor (RT) instability acting on the tangled field 
and the plasma contained within it. The tangled field acts 
like a hot, tenuous gas and is naturally buoyant; its support of 
the dense prominence plasma is likely to be unstable to flows 
that separate the gas into dense and less-dense columns. The 
observed vertical structures likely reflect the non-linear devel- 
opment of the RT instability. A detailed analysis of this insta- 



bility is compUcated by the fact that it requires internal recon- 
nection to occur within the tangled field, and it is unclear how 
rapidly such reconnection can proceed. Therefore, it is dif- 
ficult to predict the vertical velocities in prominence threads. 
Clearly, more advanced numerical simulations of the interac- 
tion of tangled fields with prominence plasma are needed. 

Finally, in Section|7]we simulated the density distribution in 
a prominence thread, using a cylindrical model for the tangled 
field and its deformation by gravitational forces. The results 
indicate that the threads have an intricate fine-scale structure. 
Multiple structures are superposed along t he LOS, consistent 
with observations of the H I Lyman lines (lOrrall & SchmahU 
[19801: iGunar etani2()07h . While our model does not take into 
account any dynamical processes, the main conclusion that 
the density distribution is highly inhomogeneous is likely to 
be valid also in the dynamical case. 

The images shown in Figures [T^ and [it were obtained at 
the Big Bear Solar Observatory, which is operated by the New 
Jersey Institute of Technology. The image shown in Figure[T]: 
was obtained at the Dutch Open Telescope, which is oper- 
ated by Utrecht University in the Observatorio del Roque de 
los Muchachos of the Instituto de Astrofisica de Canarias, La 
Palma, Spain. 



APPENDIX 
TANGLED FIELDS IN A CYLINDER 

According to the model presented in this paper, the vertical threads in hedgerow prominences are supported by tangled magnetic 
fields that pervade the dense threads and their local surroundings. In this section we construct a cylindrical model of this tangled 
field. We use a cylindrical coordinate system (r,4>,z) with r the distance from the cylinder axis (r < R where R is the cylinder 
radius), 4> the azimuthal angle, and z the height along the axis. We assume that the radial component of the field vanishes at the 
cylinder wall, Br{R,4>,z) = 0, so the field lines are confined to the interior of the cylinder. It follows that the axial magnetic flux 
must be constant along the cylinder: 

$ = / / B,(r,(f>,z)r dr d<j) = constant. (Al) 
JO Jo 

In cylindrical coordinates the force-free condition reads: 

^dB, dB^ 
r ocp oz 
dBr OB. 

We again take a to be constant (LFFF). The general solution of the above equations with the boundary condition B^iR, (j),z) = Q 
can be written as a superposition of discrete eigenmodes enumerated by an index n: 

N 

Br(r, = ^ B„Fr,„(r) sin(k„z+m„(f>+ f„), (A5) 

n=l 

N 

B^(r,(j),z) = ^ B„F^,„(r) cos(k„z +m„4> + f„), ( A6) 

n=l 

N 

B.(r,(j),z) = ^B„F,^„{r)cos(k„z + m„<j)+f„), (A7) 



n=l 



where B„ is the mode amplitude; the functions Fr^„(r), F^,„{r) and F- nif) describe the radial dependence of each mode; k„ is the 
axial wavenumber; m„ is the azimuthal wavenumber (non-negative integer); and /„ is the phase. Here n = 1 refers to a fully 
symmetric mode with ki=m\= f\=Q (also known as the Lundquist mode), while n = 2,- ■ ■ ,N refers to modes that have either a 
or z dependence. The latter modes only exist under certain conditions (see below). Inserting expressions ( IA5b . ( IA6b and (IA7l i 
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Fig. 10. — Modes of the tangled field in a cylinder (see text for details). 

into equations (IA2l i. (IA3l l and (IA4l i. we find that the radial dependencies can be expressed in terms of Bessel functions: 

Fr.nir) = -—J'^iqnr) - '^^Jmiqnr), 



F.,n(r) = JJqnr) 



qn qir 



(A8) 

(A9) 
(AlO) 



where q„ = \J o?-- k\ is the radial wavenumber, 7,„(x) is the Bessel function of order m„, and J'„j{x) is its derivative. The boundary 
condition at the cylindrical wall requkes Fr^R) = for all modes. Introducing a parameter 0„ in the range [0,7r] such that 
qn= I a I sin 6*,, and k,, = acosf?,,, we obtain the following equation for 9„: 



cos 6„J' (a sin 6',,)+ ^y„,(flsin0„) = 0, 

a sm o„ 



(All) 



where a = \a\R. The roots of this equation can be found numerically for any given values of a and m. Depending on the value of 
a, the equation may have one or more solutions: 

1. There always exists at least one solution, the axisymmetric mode (mi = 0) with 9i = tt/2. In this case qi = \a\ and ki = 0, 
so this mode is invariant with respect to translation along the z axis (Lundquist mode). 

2. There may be additional axisymmetric modes that are not invariant to translation. In the case m„ = and 6,, y^n/l, equation 
(lAl lb yields asm9„ = x\, where x/ is a root of J[pi) = (xi = 3.832, = 7.016, = 10.173, etc.). Solutions exist only 
when a > xi; for a in the range [xi,X2] there exists one such solution with d„ = arcsin(xi/fl); for a in the range [x2,ji:3] there 
exist two solution, etc. Since Jl^(a sind) is symmetric with respect to 6 = tt/2, each solution 9„ also has a second solution 
6'^ = TT-6„, but the magnetic structure of these solutions is the same, so we do not count it as a separate mode. 

3. For higher values of a there exist non-axisymmetric solutions (;«„ > 1). The first modes with m = 1 occur at a = 3.112; 
these "kink" modes are apparently more easily excited than the axisymmetric modes with ky^O. The first modes with m = 2 
occur at a = 4.708. 



For a given value of a, we systematically find all roots of equation (lAl II) . starting with m = Q and then increasing m until no more 

roots are found. Figure[TO]shows 9„ as function of a. Note that the number of modes in creas es wi th a. 

We now compute the axial flux and magnetic energy density. Inserting equations (IA7l i and (lAlOl l into equation ( lAlb . we find 
that only the Lundquist mode contributes to the axial magnetic flux: 



$ = 27rBi / Jo(\a\r)rdr=2TTRBi\a\-^Ji(a). 



where we used J[y{q„R) = for axisymmetric modes. The mean magnetic energy density is defined by 



E = 



Sir 



r dn 



(A12) 



(A13) 



where < • • • > denotes an average over (j) and z- Let nl and n2 denote two different modes, then averages of cross products such as 
< cos{k„iz+ nini(j)+ f„i) cos{k„2Z+ ni„2(t)+ f„2) > vanish. Therefore, the magnetic energy can be written as a sum over individual 



14 



VAN BALLEGOOIJEN & CRANMER 



modes, E = X]„=i E„, where 

£i = ^^ [Jl{\a\r) + Jl(\a\r)\rdr, (A14) 

^""sS^i [PrJ'-) + Finir) + Fl„{r)\rdr for«>l. (A15) 

For simplicity we consider equipartition tangled fields in which the various modes n of the LFFF have equal magnetic energy. 
This implies E„ = E /N, the same for all modes (including the Lundquist mode), which provides a relationship between the mode 
amplitudes B„. The phase angles /„ of the non-axisymmetric modes are assigned random values in the range [0,27r]. This results 
in a tangled magnetic field B(r) described by equations ( IA5b . iA6i and ( IA7I ). 
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